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ABSTRACT 

We have used the Gemini Near-infrared Integral Field Spectrograph (NIFS) to 
map the emission-line intensity distributions and ratios in the Narrow-Line Region 
(NLR) of the Seyfert galaxy NGC 4151 in the Z, J, H and K bands at a resolving 
power > 5000, covering the inner « 200 x 300 pc of the galaxy at a spatial resolution 
of ~8 pc. We present intensity distributions in 14 emission lines, which show three 
distinct behaviours. (1) Most of the ionized gas intensity distributions are extended 
to ~100pc from the nucleus along the region covered by the known biconical outflow 
(position angle PA=60/240° - NE-SW), consistent with an origin in the outflow; while 
the recombination lines show intensity profiles which decrease with distance r from the 
nucleus as / cx r"^, most of the forbidden lines present a flat intensity profile (/ cx r") 
or even increasing with distance from the nucleus towards the border of the NLR. (2) 
The H2 emission lines show completely distinct intensity distributions, which avoid the 
region of the bicone, extending from «10pc to « 60 pc from the nucleus approximately 
along the large scale bar, almost perpendicular to the bicone axis. This morphology 
supports an origin for the Il2-emitting gas in the galaxy plane. (3) The coronal lines 
show a steep intensity profile, described by / oc r~^; the emission is clearly resolved 
only in the case of [Si vii] , consistent with an origin in the inner NLR. 

Using the line-ratio maps [Fell] 1.644/1.257 and Pa/9/Br7 we obtain an average 
reddening of E(B-V) « 0.5 along the NLR and E(B-V) > 1 at the nucleus. Our line- 
ratio map [Fell] 1.257^m/[P ii] LlSQ/zm of the NLR of NGC 4151 is the first such 
map of an extragalactic source. Together with the [Fell] /Pa /3 map, these line ratios 
correlate with the radio intensity distribution, mapping the effects of shocks produced 
by the radio jet on the NLR. These shocks probably release the Fe locked in grains 
and produce an enhancement of the [Fe 11] emission at « 1" from the nucleus. At these 
regions, we obtain electron densities Ne « 4000 cm"'^ and temperatures Te ~ 15000 K 
for the [Fe ii]-emitting gas. For the H2-emitting gas we obtain much lower temperatures 
of Texc «2100K and conclude that the gas is in thermal equilibrium. The heating 
necessary to excite the molecule may be due to X-rays escaping perpendicular to the 
cone (through the nuclear torus, if there is one) or to shocks probably produced by 
the accretion flow previously observed along the large scale bar. 

The distinct intensity distributions and physical properties of the ionized and 
molecular gas, as well as their locations, the former along the outflowing cone, and 
the latter in the galaxy plane surrounding the nucleus, suggest that the H2-emitting 
gas traces the AGN feeding, while the ionized gas traces its feedback. 

Key vifords: Galaxies: active, Galaxies: nuclei. Galaxies: ISM, Galaxies: individual 
(NGC 4151) 



1 INTRODUCTION 

NGC 4151 is the nearest and apparently brightest Seyfert 
1 galaxy (or Seyfert 1.5, according to lOsterbrock fc KoskH 
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l|l976t )'). and thus harbours one of the best studied ac- 
tive R alactic nuclei (hereaf ter AGN) JCrenshaw fc Kraemerl 
I2OO7I ). As pointed out by iMundell et all ()l999[ ). from old 
optical observations, NGC4151 was believed to be a small 
spiral galaxy with major axis extending by «2.'5 along po- 
sition angle (hereafter PA) PA«130°. But this is just the 
inner part of the galaxy, comprising an oval distortion, 
or "weak fat bar" (|MundeIl fc Shonel1l999l ') . beyond which 
there is a much larger disk, weak in the optical, but bright 
in H I emission, wi th spiral arm s extending up to 6 arcmin 
from the nucleus (|Davieslll973l ). From the Hi kinematics, 
it was concluded that the galaxy has a small inclination of 
i »21°, and a major axis Vk^22° (jPedlar et al.lll992l : iDaviesI 
ll973l : lMundeU et al.|[l999l:lDas et al.ll2005l ). Its radial veloc- 
ity is cz = 997kms~^ (jPedlar et al.lll992l l. the Hubble type 
is (R')SAB(rs)ab, and we will adopt in this paper a dis- 
tance of 13.3 M pc, corresponding to a scale at the galaxy of 
65pcarcsec"^ (|Mundell et al.ll2003h . 

In radio continuum observations NGC4151 shows a 
linear structure comprising several knots elongated over 
~ 3'.' 5 along PA«77°, which is e mbedded in a dif- 
fuse emission exten ding over ~ K/.'S (|Pedlar et all 1 19931 : 
iMundell et aLlll995h Mor e recent high resolution radio im- 
ages (|Mundell et'aLll2003l l reveal a faint jet underlying the 
discrete components which seem to be shocklike features 
produced by interactions of the jet with gas clouds in the 
galaxy, as well as neutral gas absorption consistent with be- 
ing due to an obscuring torus. 

In the optical, the narrow-line region (hereafter NLR) 
of NGC 4151 has been found to have a biconical morphology 
(both in [O lll]A5007 and Ha emission lines), wit h the line-of- 
sight outside but close to the edge of the cones (|Evans et al.l 
I1993I : iHutchings et al.l[l998h . The projected opening angle 
of the cones is ~75° and the projected axis is oriented 
along a position angle PA~60/240°. Optical spectroscopy 
reveals outfl ows along the cones, with the approaching sid e 
to the SW (|Mediavilla fc Arriba j[l995l : lEvans et al.lll993h . 
The outfl ows have been modelled in a number of studies 
(e.g.JPas et a l. 2005; Crenshaw ct al. 2000; Hutclii ngs et al.l 
I1999I ). in which it is also argued that the AGN should be 
the primary ionization source of the bicone as no clear cor- 
relation with the radio jet is observed. The radio jet may be 
nevertheless associated with h igh velocity clouds observed 
in the NLR (|Winge et al.lll997l ). 

Combining X-r ay , U V an d optical spectra. 



iKraemer et"aLl (|2005l , I2OO6I ) and ICrenshaw fc Kraemerl 
1 20071 ) where able to characterize outflows estimated to be 
at only ~0.1pc from the nucleus, obtaining a high mass 
outflow rate (~0.16 M0 yr~^) which is about 10 times the 
accretion rate necessary to feed the AGN in NGC 4151. 
They sugest that this outflow originates in an accretion 

disk wind. 

In the near-infrared, iThompsonI (|l995l ) obtained spec- 
tra from 0.87/im to 2.5/xm, and concluded that, from the 
high emission- line ratios of [Fe n]/H l, the majority of the Fe 
should be in gaseous form, thus implying grain destruction 
to release a significant fraction of the iron usually tied up 
in dust. He obtained an electron temperature Te ~ lO^K 
and density A^e ~ 10* cm~^ for the region emitting [Fe 11] 
lines. Long-slit J-band observations of the kinematics of 
the [Fe ll]A1.2570/im and Pa/3 emission lines sh owed sim- 
ilar velocity structure to that observed in [Oiii] (iKnop et al.l 



Il996l ). which was later confirmed with two-dimensional In- 
tegral Field Unit (hereafter IFU) observations ([Turner et al.l 
[2002). This similarity supports the same origin for [O ill] and 
[Fe 11] emission, namely photoionization by the AGN in the 
NLR. Nevertheless, broadening of the [Fe 11] e mission rela- 
tive t o Pa/3, which is observed along the NLR (|Knop et al.l 
'199^ suggests additional processes contributing to the [Fe 11] 
emission, such as shocks from an AGN wi nd or jet. 

Despite being a well-studied galaxy l|Ulrich|[2000l ). the 
dynamics and excitation of the NLR, as well as the role of the 
radio jet, are not yet fully understood. In the present paper 
we use the Gemini Near-infrared Integral Field Spectrograph 
(NIFS) equipped with the adaptive optics module ALTAIR 
- to map the NLR gas distribution and excitation in the 
inner ~ 200 pc x 400 pc, at a spatial resolution of ~7 pc at the 
galaxy. The data cover the wavelength range 0.95-2.51 fim at 
a spectral resolving power over 5000. In a companion paper 
(|Simoes Lopes et al.ll2008h (hereafter Paper II), we use these 
data to obtain the gas kinematics and present emission-line 
channel maps of the NLR obtained by slicing the strongest 
emission- line profiles in velocity bins of 60kms~^. 

The present paper is organized as follows. In section 2 
we describe the observations and reductions. In section 3 we 
present the results, which include flux measurements of 55 
emission lines and 14 intensity distribution maps as well as 
line ratio maps. In section 4 we derive physical parameters 
for the NLR and discuss the origin of the [Fell] and H2 
emission, and in section 5 we present our conclusions. 



2 OBSERVATIONS AND REDUCTIONS 

Two-dimensional spectroscopic data were obtained on 
the Gemini North te lescope with the NIFS instrument 
(|McGregor et all 120031 ') operating with the ALTAIR adap- 
tive optics module on the nights of December 12, 13, and 
16 2006 UT. ALTAIR was used in its Natural Guide Star 
mode with optical light from the nucleus of NGC 4151 feed- 
ing the adaptive optics wave front sensor. The uncorrected 
seeing FWHM, as reported by ALTAIR, was generally in the 
range 0.6-0.9", measured in the V-band, but reached 1.2" on 
some occasions. The observations covered the standard Z, 
J, H, and K spectral bands at two-pixel resolving powers 
of 4990, 6040, 5290, and 5290, respectively. This resuhed 
in wavelength coverage of 0.94-1.16 fim, 1.14-1.36 fim, 1.49- 
1.81 ^m, and 1.99-2.42 /im, respectively. Additional spectra 
were obtained at the Kiong setting of the K grating. This 
covers the wavelength range 2.09-2.51 /im, which includes 
the H2 Q-branch. 

NIFS has a square field of view of x 3.0" x 3.0", divided 
into 29 slitlets each 0.103" wide with a spatial sampling of 
0.042" along each slitlet. The FWHM of the spatial profile 
of a star is 0'.'12 ± 0'.'02 at the H, K and Kiong bands, corre- 
sponding to «8 pc at the galaxy, while at the J and Z bands 
it is larger, 0" 16 ± 02, corresponding to ~10 pc at the galaxy. 
This is dominated by the O'.'l slitlet width across the slitlets 
and by instrumental aberrations along the slitlets. But we 
have verified an increasing strength of the uncorrected see- 
ing halo towards shorter wavelengths, which degrades the 
image quality. In order to gauge the performance of our data 
in terms of image quality, we have measured the flux in a 
0'.'2 diameter circular aperture and in a l"5 diameter cir- 
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cular aperture for each of the telluric standard stars. The 
resulting ratio between the flux in the smaller aperture to 
the one in the larger aperture is 0.19 for the Z band, 0.27 
for the J band, 0.30 for the H band, 0.47 for the K band, 
and 0.43 for the Kiong band. 

Fig-IHshows the spatial profiles of a star in the different 
bands, as compared with the spatial profiles in the galactic 
continuum in the Z, J K and Kiong bands. It can be argued 
that, as the observations were obtained using the galaxy 
nucleus as a guide star, the real PSF should be that derived 
from the spatial profile of the nuclear source, instead of the 
stellar profile. It can be seen that in the K band, the FWHM 
of the nuclear source profile is almost indistinguishable from 
that of the star, while in the Kiong band it is 0'.'04 larger 
thus 0'.'16. In the J band, the FWHM of the nuclear source 
is smaller that of the stellar profile, thus we adopt for the 
PSF the larger stellar value of 0'.'16 ± 02, while in the Z 
band the FWHM of the nuclear source is significantly larger 
than that of the star, 0'.'24. In the H band the FWHM of the 
nuclear source profile is O'.'lS ± 0'.' 02. Thus a representative 
value for the FWHM of the PSF in the J, H, K and Kio„g 
bands is 0'.'14 ± 0'.' 02, corresponding to a spatial resolution 
at the galaxy of 9± 1.3 pc. In the Z band the resolution is 
poorer, corresponding to 15 ± 1.3 pc at the galaxy. 

The instrument was set to a position angle of 345 de- 
grees to align the slitlets approxima tely perpendicular t o 
the axis of the radio jet in NGC 4151 ijMundell et al.ll2003h . 
This results in coarser spatial sampling along the jet and 
finer spatial sampling across it. The Z, J, H , and K obser- 
vations covered three adjacent NIFS fields centred on the 
NGC 4151 nucleus and offset by ± 2'.'5 along position angle 
75 degrees. The resulting field of view is 8.0" x 3.0". Only a 
single NIFS field was obtained at the Kiong grating setting 
due to the limited extent of the H2 emission. This was cen- 
tred on the nucleus at the same position angle as the other 
observations. 

Each dataset was recorded as a sequence of two 90 s 
exposures at each of the three field positions on NGC 4151 
followed by two 90 s sky exposures. The nucleus did not 
saturate in this time. The sky positions were displaced by 
« ±75" from NGC 4151 along PA = 75 degrees and dithered 
by ± 0'.'2. This sequence was repeated three times and con- 
cluded with a fourth object set. It resulted in eight object 
frames at each field position and six offset sky frames for 
each grating setting. An arc spectrum was obtained along 
with each dataset, and spectra of the nearby Hipparcos stars 
HIP 56324 (A3V) and/or HIP 61471 (AOV) were obtained 
before and/or after the NGC 4151 observations to provide 
telluric correction. The Hipparcos star observations were 
also used for fiux calibration. 

The data reduction was accomplished using tasks con- 
tained in the NIFS package, which is part of the Gemini 
IRAF package, as well as generic iraf tasks. The reduction 
procedure first applied a linearity correction and then sub- 
tracted a median-combined sky frame, multiplied by a flat- 
field frame, and cut each object frame into 29 sub-images, 
one for each NIFS slitlet. Bad pixels identified in the flat- 
field and dark frames were then removed by interpolation. 
A coordinate transformation was then applied to each two- 
dimensional sub-image to linearize the wavelength and spa- 
tial scales. These were derived from the arc exposure and 
exposures of the flatfield lamp with a Ronchi grating aligned 



perpendicular to the NIFS slitlets, respectively. The trans- 
formed two-dimensional images were then stacked into a 
three-dimensional data cube with two spatial and one spec- 
tral dimension. The data cubes for each object exposure were 
then collapsed in the spectral direction to produce a contin- 
uum image of the sky and the centroid of the NGC 4151 
nucleus was measured. The individual data cubes were then 
recentered to remove tracking drift and combined using the 
IRAF IMCOMBINE task. Each spectrum in the combined data 
cube was then corrected for telluric absorption based on the 
spectrum of the A-type Hipparcos star after intrinsic hy- 
drogen absorption had been removed by Gaussian fitting. 
The spectra were then fiux calibrated by reference to the 
telluric-corrected spectrum of the Hipparcos star, which was 
assumed to have a blackbody shape over the wavelength 
range of each near-infrared spectral band and an average 
absolute fiux density defined by the appropriate 2MASS J, 
H , or K broadband magnitude. The total signal of the Hip- 
parcos star was measured in a 1.5" diameter aperture so 
includes as much of the uncorrect adaptive-optics halo light 
from the star as practical to measure in our 3.0" x 3.0" 
field. As such, the flux calibration applies to detected light 
per spatial pixel: no correction is attempted to the total flux 
in the adaptive-optics-corrected point spread function. 

Inspection of individual spectra of telluric stars in the 
Z and J bands showed flux variations of 18% in data ob- 
tained just before and after the galaxy observations. These 
fluctuations indicate that the night of December 13, 2006, 
when these observations were made, had variable seeing, as 
there is no report of the presence of clouds during the night. 
During the previous night, when the H and K band obser- 
vations were made, and on December 16, when the Kiong 
observations were made, the seeing was stable. As a result, 
the flux-calibrated spectra of the galaxy show a flux excess 
of «17% when we compare the red end of the J band spectra 
to the blue end of the H spectra. We have thus re-calibrated 
the Z and J band fluxes dividing them by the factor 1.17, in 
order to have consistent data over all spectral bands. 

The final data cubes contain 2250 spectra per band, 
with each spectrum corresponding to a spatial coverage of 
6.6 X 2.7 pc^ at the galaxy. Although the total field covered 
in the Z, J, H , and K observations is 3" x 8", most of the 
line emission, which is the subject of the present study, is 
concentrated within the inner 3" x 5", corresponding to a 
region of dimensions 192 pc x 320 pc at the galaxy. 



3 RESULTS 

In the top panel of Fig. [1] we present a K-band image of the 
central 60"x60" of NGC 4151, obtained with the William 
Herschel Telescope (WHT), where the bar can be observed 
at PA=130°. We note that the position of the major axis of 
the galaxy (22° ) is almost perpendicular to the bar. The cen- 
tral rectangle shows the field-of-view covered by the NIFS 
observations. In the bottom panel we present an image ob- 
tained from the NIFS observations integrating the fiux of 
the [Fe 11] 1.6440/xm emission line, where we have marked 
the positions of representative spectra shown in Figs. [2] and 

m 

The spectra shown in Fig. [2] cover the Z and J bands, 
while those in Fig. |3] cover the H and K bands. From top 
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Figure 2. Sample of spectra in tlie Z and J bands, with the corresponding locations identified in the NIFS [Fen] image shown in Fig.[T] 
From top to bottom: A: spectrum from a location 1'.'4 E of the nucleus (PA=75°); B: from 0'.'7 E of the nucleus (PA=75°); C: spectrum 
from the nucleus; D: spectrum from 0'.'7 W of the nucleus (PA=255°); E: spectrum from 0''7 SE of the nucleus (PA=120°), where there 
is a maximum in the H2 emission (Fig.|6j. 



to bottom we show spectra from locations approximately at 
1'.'4 and 0'.'7 E of the nucleus (PA=75°), from the nucleus, 
from 0'.'7 W of the nucleus (PA=255°), and from C/.'? SE 
of the nucleus (PA=120°), respectively. The emission lines 
identified in the spectra are listed in Table[TJ together with 
the corresponding fluxes at the nucleus, at (/.'O SW ('.'2 to 
the right of position D in Fig.[T] where there is a peak in the 
[Fell] emission), C/.'? SE of the nucleus (position E in Fig.[T] 
where there is a peak in the H2 emission) and C/.'? E (position 
B in Fig.[IJ, for an aperture of O'.'Sl x O'.'Sl (obtained by 



binning 3 pixels along the x-axis and 7 pixels along the y- 
axis) . 

Figs.[5]and|3]show that the nuclear continuum (panel C) 
is very red, while the continuum from extra-nuclear regions 
is blue. An analysis of this continuum, and the constraints 
which can be derived for the structure producing it - pos- 
sibly adAistj'_toru^- will be presented in a forthcomig pa- 
per (|Riffel et al.ll2009l ). The nuclear spectrum shows broad 
Hi and Hel emission lines, and a deep absorption which is 
clearly observed on top of the broad profiles in the nuclear 
Hel 1.08332/im, Br7 and He 1 2.05869/im emission lines. It 
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Figure 3. Sample of H an K-band spectra, with tlie corresponding locations identified in the NIFS [Fell] image shown in Fig.[T] and 
specified in the caption of Fig|2] 



can also be observed that there is no H2 emission at the 
nucleus and that the coronal lines of [Ca vill] and [Si vil] 
seem to be extended. The [P 11] 1.1886/^m emission line is 
also extended and can be observed up to l'.'4 E of the nu- 
cle us. This line was fi rst reported in an extragalactic source 
by lOliva et al] l|200ll ). and since then observed in the nu- 
cleaj spectrum of le ss than 15 galaxies, including NGC4151 
l|Riffel et all 1200681 ). Combined with [Fe n]1.2570/im, the 
[Pii] line is a powerf ul diagnostic of the origin of the [Fe 11] 
emiss ion in galaxies l|01iva et al.l I2OO1I : Ijackson fc Beswickl 
I2OO7I ). In the present paper we provide the first 2D map of 
an extragalactic source in this emission line. 



The profiles of most of the emission lines vary accord- 
ing to the location in the NLR, and an illustration of this 
variation is shown in Fig.|4j where we present a sequence 
of selected emission-line profiles, from different locations 
along the bicone (PA=60°). A vertical dashe d line shows 
the a dopted systemic velocity, of 997kms~^ jPedlar et al.l 
Il992h . The same fiux scale is kept for each emission line 
in Fig|3] so that the fiux variation can also be observed 
as a function of distance from the nucleus. The scale nev- 
ertheless varies for the different emission lines. It can be 
observed that the center of the emission lines drifts from 
blueshifts observed to the SW of the nucleus to redshifts 
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Figure 4. Sample of profiles of the [S III] 0.9533/im, He 1 1.0833^m, [Fe ll] 1.2570 /im and Br72.1661^m emission lines, from spectra 
extracted along the bicone axis (PA=60°), at the locations relative to the nucleus indicated in the figure. The same flux scale is kept at 
all locations for each emission line, but varies from line to line. 
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Table 1. Emission line fluxes within 0'.'3 X 0'.'3 apertures (10 erg cm s ^) 



A (vac.) ID Nucleus [Fc ll] peak (O" 9 SW) H2 peak (O" 7 SE) Pos. B(0"7E) 
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observed to the NE of the nucleus. In some locations - e.g., 
at O'.'S SW of the nucleus in Fig|4]- many emission lines are 
clearly double-peaked, while in other locations the profiles 
are asymmetric suggesting the presence also of two com- 
ponents, although unresolved. The blueshifted absorption 
in the Hel 1.08332/im emission-line profile is observed at 
— 468kms~^ relative to the systemic velocity (more details 



will be given in Paper II). A similar absorption (with simi- 
lar blueshift) is observed in Br7 and in the Hel2.0587^m 
lines. The absorption in both He I lines reach below the 
interpolated continuum indicating that the absorber cov- 
ers much of the continuum source and the BLR clouds. 
To our knowledge, these absorption features have not been 
seen in previous near-infrared spectra, such as the one ob- 
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WHT (K-band) 
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Figure 1. Top: K-band image of the central 60"x60"of 
NGC4151 obtained with WHT. The image has been rotated to 
the same orientation of the NIFS frames. The continuous hne 
shows the orientation of the major axis of the galaxy, the dashed 
line shows the orientation of the bicone and the dot-dashed line 
shows the orientation of the bar (visible in the image). The rect- 
angle shows the region covered by the NIFS observations. Bot- 
tom: image obtained from the NIFS frames integrating in the 
[Fell] A 1.644/.im emission line; the letters identify locations corre- 
sponding to the spectra shown in Figs .[2] and |3] 



tai ned bv [Osterbrock. Shaw fc Veilleuxl (jl990l ) in 1988 or 
bv iThompsonI (|l995h in 1993. More recently, iRiffel et al.l 
l|2006a ) also fail to detect this absorption. Nevertheless, 
absorptions of similar widths and blueshifts have b een ob- 
served in the UV and optical. iKraemer et al.l (|200ll) report 
absorptions at a velocity — 490kms~^ with respect to sys- 
temic in HST UV spectra of the NLR of NGC4151, which 
had been previously identified as d ue to component s calle d 
D an d i? in a previou s study by IWevmann et al.l |l993), 
while iHutchings et al.l (2002) report an absorption due to 
HelA3889A at -460kms"\ The absorptions we see in the 
near-IR are probably related to these UV and optical ab- 
sorptions seen in HST spectra. Its detection in our spectra is 
aided by the high spatial resolution of the NIFS data, which 
minimizes NLR contamination of the nuclear spectrum. 



3.1 Emission-line intensity distributions 

2D maps of the emission-line intensities have been obtained 
by integrating the flux under the line proflles, after subtrac- 
tion of the contribution of the underlying continuum, deter- 
mined as the average between two spectral windows adja- 
cent to the emission lines. In the case of the H and He lines, 
which have broad components at the nucleus, the adjacent 
continua faU on top of the broad lines and thus the resulting 
flux is essentially from the narrow component. Nevertheless, 
as the broad-line proflles are not symmetric (e.g. Pa/3 proflle 
in the middle panel of Fig.[2)| the subtraction is not always 
perfect and the nuclear fluxes in these lines may have some 
residual broad-line flux. 

The intensity maps for the main extended 
emission lines from the Z and J bands, namely 
[S III] O.QSaS^tm, [Sviii] 0.9915 /xm. He ii 1.0126Atm, 
[Sii] 1.029,1.032,1.034,1.037Aim, He: 1.0833, [P ii] 1.1886^m, 
[Six] 1.2523^im and [Feii] 1.2570 ^im are shown in Fig.[5l 
Along the y-axis, we show the whole extent of the frames 
(3"), while along the x-axis we show only the region 
with measurable flux in at least one of the emission lines, 
comprising the inner 4'.' 6. We have overplotted on the He: 
flux map the contours (in blue) of a radio A21 cm image 
from iMundell et all (|l995l '). obtained with MERLIN. These 
contours have been aligned to the intensity maps under 
the assumption that component C4 of the radio image is 
aligned with the peak of our K-band continuum map. Radio 
component C4 is claimed to contain the active nucleus. 
We assume that the K-band continuum peak also contains 
the nucleus. The continuous line in the top left panel 
shows th e orientation of th e major axis of the galaxy at 
PA=22° (jPedlar et al.lll992l ). while the dashed line shows 
the orientation of the [O ill] bicone at PA=60/240° and 
the dot-da shed line shows the ori entation of the bar, at 
PA=130° l|Mundell fc Shond 1 19991 ). The [P 11] flux map 
shown in the second panel (from top to bottom) to the 
right of Fig. [5]), is the flrst 2D map in this line ever obtained 
for an extragalactic source (as pointed out in the previous 
section). The resolution achieved in the images can be 
judged by a comparison with telluric star images, show as 
insets in Fig.[S] the image of a star in the Z band is shown 
in the [S ill] panel, while that in the J band is shown in the 
[P 11] panel, with peak intensity normalized to that of the 
galaxy in the emission line of the corresponding panel. 

In Fig.[6]we show the intensity maps for additional emis- 
sion lines observed in the J, H and K bands, namely Pa /3, 
[Fe 11] 1.6440 /xm, H2 2.1218 ^m, Br 7, [Ca vill] 2.3220 and 
[Si vil] 2.4833 pim. We have overplotted (in blue) on the [Fe 11] 
map the contours of the [O ill] emission-line map. 

From Figs. [5] and [S] it can be observed that the light 
distribution varies for different emission lines. For the coro- 
nal lines [S vill] , [S ix] and [Ca vill] , the light distribution is 
compact and not clearly resolved, while the coronal [Si vil] 
emission distribution is more extended, being aligned with 
the bicone to the SW and with the radio jet to the NE. 

Except for H2, the other emission lines show much 
stronger extended emission, reaching at least 4" along the 
bicone and 2" in the perpendicular direction, correspond- 
ing to projected distances at the galaxy of 256 pc and 
128 pc, respectively. The intensity distributions are some- 
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Figure 5. Intensity maps of emission lines from the Z and J bands, obtained by integrating tlie flux under the line proflles after 
subtraction of the continuum. The green contours correspond to 3 times the background noise (Sit). The blue contours are from the radio 
MERLIN image of lMundell et alj lll995l') . The dashed line on the top left panel shows the orientation of the bicone, the continuous line 
shows the galaxy major axis, and the dot-dashed line shows the orientation of the large scale bar. The insets show images of telluric stars 
in the Z band (in the [S in] panel) and J band (in the [P ll] panel) with peak intensity normalized to that of the galaxy in the emission 
line of the corresponding panel. Flux units are 10~^^ ergcm"^ s~^ spaxel"^. 
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Figure 6. Intensity maps of emission lines fro m the J, H and K ban ds, represented as in Fig. [5] Blue contours overplotted on the 
[Fell] A 1.644 fim map are from the [Olll] image of lHutchings et"aL I lll998l '). The insets show images of telluric stars in the H band (in the 
[Fen] panel), K band (in the H2 panel) and in the Kjo^g band (in the [Sivil] panel), normalized to the peak intensity of the galaxy in 
the emission line of the corresponding panel. 



what brighter and more extended to the SW (the near side 
of the bicone) than to the NE (the far side). 

The [O III] image, whose contours are overplotted on 
the [Fell] intensity distribution in Fig. [6] has a "knotty" ap- 
pearance, similar to what is observed in our [Fell] map, with 
some knots coinciding in both maps. Nevertheless, a detailed 
comparison shows that the decrease in emission just outside 
the nucleus and increase again in two opposite regions at 
about 1" from the nucleus along the bicone observed in the 
[Fe 11] maps is not observed in the [O ill] map, which shows 
emission in "curv ed strands" all the way along the bicone 
iHutchings et al.| [r998). 

A completely different light distribution is observed in 
the H2 flux map (bottom left panel of Fig.|6]): instead of 



being elongated along the bicone axis, it is elongated al- 
most perpendicular to it. There is almost no H2 emission 
within «0'.'3 from the nucleus, its flux being distributed in 
two structures resembling double arcs to the NW and SE, ex- 
tending from ~0'.'3 to ~l"from the nucleus along the minor 
axis of the galaxy (PA=112/292°), which is also approxi- 
mately the orientation of the bar. 

In order to look for a possible relation between the 
emission-line intensity distributions and the radio emission, 
we have ove rplotted the radio con tours of the MERLIN ra- 
dio image of lMundell et al.1 l|l995h on the He l flux map (top 
right panel of Fig.[S]). While the radio intensity distribution 
is oriented horizontally in the Figure, along PA=75°, most 
intensity distributions are oriented instead along the [O ill] 
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Figure 7. Line rat io maps, where the green contours overplotted on the [Fell]/Pa/3 map are from t he radio MERLIN ima ge of 
iMundell et al] lll995l) . while the ones overplotted on the [Fe ll]l. 644/1. 257 map are from the [Olll] image of lHutchings et aL 
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bicone at P.A. =60/240° . Thus there seems not to be a strong 
association between the radio jet and the hne emission. On 
the other hand, both the radio and hne emission are stronger 
to the SW than to the NE, and a closer inspection shows 
that the radio component C3 (the radio knot just W of the 
nucleus) seems to align well with a region of enhanced emis- 
sion 0'.'4 W from the nucleus observed in the [Fe 11] 1.6440/im 
map. 

3.2 Line ratio maps 

We have used the intensity maps to construct the line 
ratio maps [Sn]/[Sni], [Sni]/Pa/3, [Feii] 1.257/Pa/3, 
[Fell] 1.257/ [P II], H2 2.122/Br7, [Fe 11] 1.644/1.533, 
[Fell] 1.644/1.257 and Br7/Pa/3, which are shown in 
Fig-El As pointed out above, at the nucleus (and within 
!«0'.'2 from it), there may be some contribution from the 
broad-line component in Pa/3 and Br7, thus the nuclear 
ratios involving these two lines may be affected by this 
component. Many of these lines ratios are indicators of the 
gas excitation, thus the nuclear values should be used with 
caution. 

The [S III] /Pa/3 line- ratio map has values of « 4 at the 
nucleus, increasing outwards to ~ 12 to the SW and r; 8 to 
the NE. There is thus a systematic difference in the values 
of this line ratio between the two sides of the bicone: higher 
ratios in the near side (SW) and lower ratios in the far side 
(NE). The [Sll]/[Slll] ratio map has the lowest values <0.1 
at the nucleus, increasing outwards up to 0.2. 

We have overplotted the contours of the MERLIN radio 
image on the line ratio [Fe 11] /Pa/3 map in order to look for a 
possible relation between the [Fe 11] excitation and the radio 
emission. There is indeed a relation: both the [Fe ll]/Pa/3 and 
[Fell]/[Pll] (Fig.[7| line ratios increase outwards, reaching 
maximum values at «l" SW, the location where the radio jet 
shows a "flaring" in its distribution, and to the opposite side 
at ?»1" NE, the location of a radio hotspot. The [Fell]/Pa/3 
ratio increases from values < 1 close to the nucleus up to « 3 
at the locations of the radio SW flare and NE hotspot, while 
[Feii]/[Pii] increases from values ~ 4 to « 8 at these same 
locations. The line ratio [Fell] 1.644/1.533, which is sensitive 
to the emitting gas density, shows values < 4 within the 
inner 0'.'5, increasing to ~ 8 — 10 in the outermost regions. 

The H2 /Br 7 line ratio distribution is a result of the H2 
intensity distribution, which shows low or zero fluxes at and 
close to the nucleus and higher values perpendicular to the 
bicone. The consequence is an increase in the line ratio from 
values ^1 around the nucleus to >3 towards the borders of 
the H2 emitting regions at ~ l"to the NW and SE of the 
nucleus. As in the case of Pa/3, within O'.'S from the nucleus, 
the line ratios may be affected by the broad component of 
the Br7 line, and thus should not be used as indicator of the 
gas excitation. 

The line ratios [Fe 11] 1.644/1.257 and Brj/Pa.f3 are 
shown in the bottom panels of Fig. [7] These line ratios can 
be, in principle, used for estimates of the reddening along 
the NLR, although the wavelength baseline is small, mainly 
for the [Fe 11] ratio. The highest values are observed at the 
nucleus, for both ratios. In the [Fell] 1.644/1.257 ratio map 
there is a hint of a structure resembling the curved strands 
observed in the [O ill] image (whose contours are overplot- 
ted on this map) in which the highest ratios seem to be ob- 



served in the regions of lowest [O ill] fluxes. Data with better 
spatial resolution and higher signal-to-noise ratio would be 
necessary in order to confirm this result. 



4 DISCUSSION 

4.1 Ionized gas distribution 

The intensity distributions in the emission lines of the 
ionized gas (with the exception of the [Svill], [Six], 
[Ca vill] and [S vil] coronal lines) resemble that of the 
optical [Olll] emission line, suggesting a similar origin, 
namely emission from ionized gas outflowing along the 
walls of a hollow bicone centere d at the nucleus and ori- 
ented alo ng PA=60/240° (Hutchings et all 1 19991 : [Pas et all 
2005; Cre nshaw et al boootT Nevertheless, as pointed out by 
iKraemer et al.l ( 20081 ), these intensity distributions do not 
clearly delineate a bicone, showing also extended emission 
perpendicular to the axis of the bicone in the vicinity of the 
nucleus. In other words, there is line emission beyond the 
presumed walls of the bicone close to the nucleus, indicating 
some escape of radiation in the perpendicular direction. 

The fact that the light distributions are more extended 
and somewhat brighter to the SW than to the NE can be 
understood as due to the SW cone being directed toward us 
and we are looking inside the cone, where we are observing 
the gas most exposed to the nuclear radiation fleld. To the 
NE we are observing the emission through both the wall of 
the far cone and the galactic plane (see discusssion below 
and in Paper II). 

Radial proflles of the NLR emission have been obtained 
by averaging the line flux within conical regions having an 
opening angle of 40° centered on the nucleus and oriented 
along the bicone axis at PA — 60° and 240°. These ra- 
dial proflles are shown in Fig. |8l Profiles of [0 11] 3727A 
and [O III] 5OO7A are also plotted from data presented by 
iKraemer et al.l |20o3). All profiles decline steeply beyond 
the edge of the bright inner NLR at radii between 1" and 
2". Among the near-IR profiles, the strongest features are 
those of [S III] 0.9533 /im (green crosses) and He I 1.083 /im 
(black open circles). 

Three distinct behaviours can be observed for the radial 
profiles in Fig[S] The first is observed for [S ill] 0.9533 ^m 
(green crosses) and the recombination lines of He 1 1.083 /im 
(black open circles), Pa/3 (black crosses) and Helll.013/^m 
(dark blue crosses). Their intensity profiles decrease mono- 
tonically with distance from the nucleus, showing a depen- 
dence on the radial distance of « 7 oc , as illustrated 
by the black dashed lines in Fig.[S] In order to understand 
the behavior of the lines, it would be necessary to con- 
struct models to try to constrain the physical parameters 
which lead to line emission, what is beyond the scope of the 
present study. A pre liminary comparison of the radial pro - 
files with models by Groves. Dopita. fc Sutherland! (|2004ah 
an d[G roves. Dopita! fc Sutherland! (|2004bD suggest that pro- 
files such as the above are reproduced by models including 
dusty radiation-pressure dominated clouds. Dust-free mod- 
els for the NLR result in radial proflles for the Pa /3, He I and 
He II emission lines which are steeper than those observed. 

A distinct behaviour is observed in the [Fe 11] radial 
proflle (red crosses in Fig.|8}, which is enhanced at radii 
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around ~ 1" due to the emission clumps that are apparent 
in Fig. [6] This could be explained by excess gas-phase abun- 
dance caused by shocks produced by the radio jet that has 
destroyed dust grains and released Fe. Support for this inter- 
pretation is given by the correlation between by the [Fe ll] 
/[P ll] ratio and the radio jet. In Paper II, we also find a spa- 
tial correlation between the radio jet and the [Fe ll] kinemat- 
ics. The behavior of [Fe ll] emission can also be partly due to 
the fact that [Fe ll] is produced in partially ionized regions 
beyond the main hydrogen ionization front in NLR clouds. 
Progressive absorption of FUV photons near the hydrogen 
ionization edge by absorbing clouds located between the cen- 
tral black hole and the NLR hardens the photoionizing spec- 
trum. It is possible that the NLR has become optically thick 
to hydrogen-ionizing FUV radiation in its outer parts while 
it remains optically thin to X-rays from the AGN. In these 
circumstances, the outer NLR clouds could develop exten- 
sive partially ionized regions that emit relatively more [Fe ll] 
than P /?. 

[O ll], [O III] and [Sll] (black asterisks, black filled cir- 
cles and red circles, respectively) show a similar behaviour 
to that of [Fe ll] , which can be approximately described as 
/ cx r°. 

A third behaviour is observed in the coronal lines, which 
we now discuss. 

4.2 Coronal gas distribution 

Fig. |8] shows that the coronal lines (light blue and magenta 
symbols) all decline more steeply with radius than the re- 
combination lines. In fact, they are well-represented by an 
/ (X intensity decline, similar to that observed for the 
star radial profiles (shown as dashed blue and red lines in 
Fig.[SJ. The [Cavill] feature (light blue circles in Fig. 
declines even more steeply at small radii. 

In order to verify if the coronal emission is resolved 
and quantify its extent, we have constructed azimuthal av- 
erages within circular radial annuli of the light distribution 
in each coronal line and have normalized the profile to unity 
at the peak. The resulting spatial profiles are shown in Fig. [5] 
together with those in the galaxy continuum and telluric 
standard stars. In the top panel we compare the profiles 
of the [S vill] A0.991 /im and [S ix] A 1.252/im coronal lines 
with those in the continuum and telluric star in the Z and 
J bands, respectively. In the bottom panel, we compare the 
profiles of the [Ca vill] A2.322/im and [Si vil] A2.483/im coro- 
nal lines with those in the continuum and telluric star in 
the K and K;o„g bands, respectively. As discussed in Sec. [51 
if we assume that the PSF is given by the profile in the 
galaxy continuum, instead of the stellar profiles (except for 
the J-band, where we adopt the broader stellar profile), one 
concludes that none of the coronal lines is strongly resolved, 
except for [Si vil], which is somewhat extended along the bi- 
cone axis, what can be seen already in Fig.jS] [Cavill] may 
be slightly more extended than the nuclear continuum, but 
higher spatial resolution and better sampled data, possibly 
with better signal-to-noise ratio, would be required to con- 
firm this given that the NIFS spaxels have a size of O'.'l in 
one direction. 

The fact that the light distributions in [S ix] and pos- 
sibly [Cavill] coronal lines are extended, support an origin 
for these coronal lines in the inner part of the NLR or in 
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Figure 8. Radial emission-line profiles in cones of 40° opening 
angle at PA = 60° (top) and 240° {bottom) along the bicone axis; 
[Fell] 1.257 tim (red crosses), [Sll] 1.029/32/33 fim {red open cir- 
cles), Hi P/3 1.282 /xm (black crosses), [Slll] 0.9533 /^m {green 
crosses). He I 1.083 fim {black open circles), Hell 1.013 /im {dark 
blue crosses), [On] 0.3727^m {black asterisks), [O III] 0.5007/im 
{black filled circles), [Sivil] 2.483 fim {light blue crosses), [Cavill] 
2.322 /im {light blue circles), [Svill] 0.991 {magenta crosses), 
[Six] 1.2523 ^tm {magenta open circles). Telluric star profiles in 
the J and K bands are shown in blue and red dashed lines, re- 
spectively. Black dashed lines illustrate the slopes corresponding 
to diff'erent radial dependences for the intensity. 
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Figure 9. Comparison of the spatial profiles of a star in the J 
and K-bands (dashed lines) with those of the galaxy nucleus in 
the J and K continuum and in the coronal lines [Six], [Svill], 
[Caviii] and [Sivii]. 



Table 2. Ionization Potentials 



Line 


IP (eV) 


[Fell] 


7.9 


[Sll] 


10.4 


[Pll] 


10.5 


[Siii] 


23.3 


[Ca viii] 


127.2 


[Sivii] 


205.3 


[S viii] 


280.9 


[Six] 


328.8 



the transition region between the BLR and the NLR, as 
sugge sted by previous authors (e.g. iRodriguez-Ardila et al.l 
I2OO6I ). Ionization potentials for the parent species of each 
coronal line emitter are listed in Table [2] These gener- 
ally support the formation of the observed stratified coro- 
nal emission region through photoionization by the central 
AGN, with higher ionization potential species located closer 
to the ionizing source. With photon energies of 127.2 eV 
being required to create Cavill, the smaller extent of the 
[Cavill] emission region when compared to that of the 
[SiVIl], which has higher ionization potential is unexpected. 
This fact can be attributed to depletion of calcium onto dust 
grains, which can alter the gas-phase calcium abundanc e by 
large factors (e.g. iGroves. Dopita. fc Sutherland|[2004al lbl). 



panel of Fig. [6]) and is more extended than other emission 
lines along the minor axis of the galaxy, which coincides ap- 
proximately with the orientation of the bar. The H2 emission 
covers a region fa 20-60 pc in radial extent from the nucleus. 
One possible explanation for this morphology is that the nu- 
clear molecular gas is located predominantly in the plane of 
the galaxy, but is dissociated by the AGN radiation field 
within the bicone. It is known that the bicone intercepts 
part of the galaxy plane Q. Perpendicular to the bicone, the 
H2 must be shielded from the strongest dissociating radia- 
tion, probably by a dusty torus and/or by the walls of the 
bi-conical outflow l|Kraemer et al.ir2008l '). 

A previous H2A2.1218/im image of the nuclear r egion 
of NGC4151 was otained by iFernandez et all l|l999l ) who 
have found a similar intensity distribution to ours, in the 
form of a partial ring surrounding the nucleus, which led 
them to propose that the H2 emission originated in the 
outer part of a molecular torus. They suggested that the 
H2 emitting gas could be rotating in a plane perpendicular 
to the radio axis. We note, however, that th e orientation 
of the partial ring shown in Fernan dez et al.l ([l999) is ro- 
tated by 90° and is mirrored relative to our H2 intensity 
distribution. Otherwise, the intensity distributions are con- 
sistent with each other considering t he di fference in spatial 
resolution (0'.'6 in [Fernandez et al. I (|l999t ) and O'.'ll here). 
As our datacube gives intensity distributions for the ionized 
gas in agreement with that for the well known bicone, and 
the H2 intensity distribution is obtained from the same dat- 
acube, we c onclude that the orient ation of the H2 ring is 
mistaken in 'Fernandez et all ^199^ ). Our observations also 
do not support the suggestion by these authors that the H2 
gas is rotating in a plane perpendicular to the radio jet. In 
Paper II we present the H2 kinematics which shows very lit- 
tle rotation consistent with a gas distribution in the galaxy 
plane along the minor axis. Nevertheless, our observations 
are consistent with the interpretation that the H[2 emitting 
gas may be tracing the gas reservoir which feeds the SMBH. 
Results supporting this idea are the inflows measured by 
[Mundell & Shone (1999) in radio observations of Hi along 
the large scale bar (at PA=130°). These inflows direct gas 
towards the inner part of the bar in the nuclear region and 
their orientations are approximately perpendicular to the 
arc-shaped structures delineated by the H2 light distribu- 
tion (see Fig.[6ll. Thus one possibility is that the radio ob- 
servations are tracing H I inflows towards the nuclear region 
leading to the build up of the molecular gas reservoir which 
we observe in the H2 intensity distribution. 



4.4 Extinction 

Previous stud ies report low redden i ng to wards the NLR 
of NGC4151. ICrenshaw fc Kraemerl (|2005l ) claim only low 
Galactic extinction of E{B - V) = 0.02 ± 0.04 mag in 
STIS spectra c e ntere d on the "nuclear emission-line knot". 
iKraemer et al.l l|2000t ) report that at other locations in the 
NLR the optical reddening of the emission lines ranges from 
E{B-V) ^ 0.0 to 0.4. 



4.3 Molecular gas distribution 

The H2 intensity distribution is totally different from those 
of the other emission lines. It avoids the bicone (bottom left 



^ See Paper II and iDas et al] (e.g. l2005h for a comprehensive 
discussion of the kinematics and geometry of the bicone and the 
galaxy. 
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Other studies (e.g. [Alexander et al.lll999l ) report red- 
dening estimates which range from almost negl igible E(B — 
V) = 0.04-0.05) (e.g.JKriss et al. 1995 : Pcnston et al.lll98ll l 
to the considerable E{B ~ V) = 0.13 (|Malkanlll98i r From 
a vari ety of methods (an d data), including near-IR emission 
lines. IWard et all (Il987l) get E{B - V) ~ 0.23 mag, while 
iRieke fc Lebofskvl_(ll98l1') o btain 0.5 < Av < 0.8. 

Mundell et al.1 (|l995h measure H I absorption to 
the radio component C 4 of 3.9 x 10^^ cm^'^. From 
iBohUn. Savage, fc Drake! (119781 ) < N{HI + H2)/E{B - 
V) >= 5.8 X lO^'^ cm~^ mag~^, so this Hi column corre- 
sponds to E{B ~ V) = 0.7 {Av 2.1 mag) to the nuclear 
radio source, if there is no intervening H2. 

We can use the line ratio map [Fell] A1.2570/A1.6440 to 
estimate the reddening to the NLR of NGC 4151. These lines 
arise from the same upper level so the intrinsic value of the 
line ratio s hould be 1.36, according to t he transition prob- 
abilities of iNussbaumer fc StorevI ll 19881) a n d as confirmed 
from observations by Bautista fc PradhanI (|l998l ). Smaller 
line ratios indicate the presence of reddening and can be 
used to estimate its value through the relation: 



I I I I I I I I I I I I I I 



NE (SE) arcsec (NW) SW 



-1 1 

NE (SE) arcsec (NW) SW 



E{B ~V)= 8.14 X 



1.36 

^1.2570 
^1.6440 



(1) 



of 



obtained using the reddening law 

ICardeUi. Clayton fc Matli3 (igsi). 

Typical values of the reddening uncertainty can be ob- 
tained from the data in Table[T] The measured flux ra- 
tios are 1.21 ± 0.06 and 1.25 ± 0.12 for the SW [Fell] 
peak and Position B, respectively. These ratios equate to 
E{B -V)= 0.41 ± 0.14 and 0.30 ± 0.28, from which we 
estimate random uncertainties in E{B — V) derived by this 
method to be ~ ±0.2 mag along the radio axis. For the nu- 
cleus and other regions where the [Fe 11] emission is weaker, 
the uncertainties are twice as large. 

In order to investigate further the reddening variation 
across the NLR, we have extracted smoothed ID profiles 
from the [Fell]A1.644 and 1.257/^m intensity maps along a 
pseudo slit of width 0'.'3, oriented along the bicone. The 
fluxes were obtained at each O'.'l along the pseudo slit as 
the average of the fluxes of the pixels included within the 
O'.'S width of the slit. The ID profiles were then smoothed 
further by replacing the flux at each position by the average 
of its flux and those of the two adjacent positions along the 
slit. We have then constructed the ratio between the two ID 
profiles along the bicone and obtained the E{B — V) from 
the expression above. The result is shown in the right upper 
panel of Fig. 1101 showing a maximum E{B — V^) = 1.4 at 
the nucleus, which decreases abruptly to values in the range 
0.3 < E{B - V') < 0.6 beyond 0'.'2 from the nucleus, with 
an average value of E{B — V)!^ 0.45. 

In order to check the above result, we have next used 
the Pa /9/Br 7 ratio to estimate the reddening along the same 
pseudo slit, using the relation: 



E(B - V) = 4.74 X 



(2) 



where we have used again the reddenmg law 01 
ICardeUi. Clavton fc Mathid /l989^ and adopted the intrinsic 
Fpis/Fbt-i ratio of 5.88 corresponding to case B recombina- 



Figure 10. Left: spatial line ratio profiles along the bicone axis 
(PA= 60°) obtained from smoothed ID profiles extracted within a 
pseudo slit of width O'.'S. Filled circles show values not corrected 
for reddening, while the open circles show values corrected for 
reddening with E(B-V) obtained from the average of the two E(B- 
V) values from the right panels for each location. Right: E(B-V) 
values obtained from the line ratio profiles of [Fe 11] Al. 257/1. 644 
(top) and Pa/3/Br7 (bottom). 



tion (|Osterbrocklll98^ ). Due to the somewhat longer wave- 
length baseline, the uncertainties in the resulting E{B — V), 
of ±0.15, are smaller than that obtained using the [Fell] 
ratio. The results are shown in the bottom right panel 
of Fig.[TOl Within the uncertainties, the E{B — V) values 
and spatial variation are similar to those obtained using 
the [Fell] ratio: the highest values are observed at the nu- 
cleus, falling abruptly outwards to an average value around 
E{B-V) ^0.^. 

We note that the nuclear E( B-V) values ar e in ag ree- 
ment with the one obtained by iMundell et al.l (|l995l ) for 
the nuclear source, while E(B-V)»i0. 5 is consistent with th e 
highest E(B-V) values obtained by iKraemer et al.l l|2000l ), 
considering also the fact that the optical depth reached in 
the near-IR is larger than those obtained by optical obser- 
vations. Finally, we pointed out at the end of Section[33] 
that there may be a relation between the [Fe 11] Al. 257/1. 644 
map and the [O ill] intensity distribution. In terms of red- 
dening distribution, this relation would imply that the re- 
gions with largest reddening would coincide with the gaps 
in the intensity distribution, suggesting that the appearance 
of the [O III] intensity distribution could be, at least in part, 
due to reddening. As E(B-V)~0.5 has been obtained from 
smoothed profiles, this value is consistent with higher ex- 
tinction in narrow structures not resolved by these smoothed 
profiles. For example, if the average E(B-V)«0.5 is due to 
the smoothing of narrow strands having E(B-V) = 1, with 
E(B-V)=0 between the strands, the [Olll] flux will be re- 
duced by a factor of fti25 where E(B-V) = 1, what could result 
in the stranded appearance observed in the [O lll] intensity 
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distribution. Data with better spatial resolution and higher 
signal-to-noise ratio would be necessary in order to test this 
possibility. 

4.5 Gas excitation 

The line ratio [S ill] /Pa/? (top right panel of Fig. 
O can be used as a tracer of the excitation of 
the NLR gas. We have obtained typical values for 
this line ratio in active and starburs t galaxi es from 

fc ChaUid 
case of 
we have 



the works 
i 19951) and 



of IStorchi-Berpjinann. K inney 



Riffel et al 



2006al\ . 



In the 

(Eiii), 



IStorchi-Bergmann. Kinney fc Challid 
used their observed [S ill] /Ha ratio values (not corrected 
for reddening) to estimate the [Slll]/Pa/3 ratios (assuming 
a typical reddenin g of E(B-V)=0.5), while in the case of 
iRiffel et all l|2006ah we have used their measured values for 
[S III] and Pa/3 fluxes, also not corrected for reddening. In 
order to avoid the contribution of the broad components 
to the emission lines, we have collected only the values for 
Seyfert 2 galaxies as typical of the NLR, as well as those 
for starburst galaxies, to use as a comparison. The range 
of line ratios obtained for a total sample of 14 Seyfert 2 
galaxies is 6 ^ [S lll] ^ 14, with most values clustering 
around [Slll]~8.5, while for a sample of 10 starburst 
galaxies the observed range is 1 ^ [Sm] 5; 3.5.. These 
values can be compared to those in Fig[7l In order to make 
such a comparison easier, we present, in the bottom panel 
of Fig llOl line ratio profiles along the bicone axis obtained 
from smoothed ID profiles extracted along a pseudo slit of 
width C/.'S, as described in the previous section. 

In the ID profiles of fig.[TOl the observed [S lll]/Pa/3 line 
ratios (filled circles) beyond 0'.'5 from the nucleus are in the 
range 7-10, thus similar to the values for NLR of Seyfert 2 
galaxies, with the values to the NE being 20-30% smaller 
than those to the SW. Within O'.'S from the nucleus the line 
ratio decreases down to half the values beyond this region. 
This decrease is partially due to reddening, as can be ob- 
served in the reddening corrected line ratio profile shown 
as open circles in the bottom left panel of Fig. 1101 The cor- 
rected values were obtained by us i ng ag ain the reddening 
law of ICardeUi. Clayton fc Mathi^ (| 19891 ) and the E(B-V) 
for each location as the average of the two values shown 
in the right panels of Fig. 1101 The values of the line ratio 
within 0'.'5 from the nucleus increase from 4-6 to 7-8, now 
within the typical range observed for AGN. Nevertheless, the 
correction for reddening does not eliminate completely the 
decrease observed at the nucleus, and another effect may be 
present, possibly a residual contribution of the broad com- 
ponent of H/3. In principle, the 20-30% smaller [Slll]/Pa/3 
values to the NE relative to the SW could also be due to an 
excess reddening of E(B-V)~0.5. Nevertheless, the E(B-V) 
profiles in the right panels of Fig. [10] do not show a system- 
atic difference in reddening between the SW and NE of this 
order, and we thus attribute the lower ratios to NE to lower 
excitation. 

The [S ii]/[S III] ratio (top left panel of Fig. [T]! can also 
be used as a tracer of the excitation, although this ratio is 
noisier than the [Siii]/Pa/3 ratio due to the fact that the 
[S ll] flux map is a collection of the fluxes in four faint emis- 
sion lines. We have also constructed a ID line-ratio spa- 
tial profile for [S ll]/[S ill], which is shown in the top left 



panel of Fig. 1101 Beyond 0'.'5 from the nucleus, this ratio is 
higher to the NE, which indicates lower excitation relative 
to the SW side of the bicone, in agreement with the be- 
havior of the [S iii]/Pa/3 ratio. Within 0'.'5 from the nucleus, 
the [S ii]/[S III] ratio decreases, what cannot be attributed to 
reddening, as an excess reddening should increase the value 
of this line ratio and not decrease. Although this decrease 
could be due to a higher excitation in this region, this is 
not consistent with the behavior of the [S lll]/Pa/3 line ratio. 
Our favored explanation is as follows. In the next section 
we use the observed [Fe ll] line ratios to derive the electronic 
density along the NLR. Within « 0'.'5 from the nucleus, the 
average density (see Fig. [12] is higher than the critical den- 
sity of [S ll] (2.5x lO^cm"^), resulting in a decrease in the 
intensity of the [S ll] lines close to the nucleus. 

How can one understand the lower excitation to the 
NE as compared to t he SW side of the cone? According to 
previous models fe.g. j Hutchings et al ] | 19981 : ICrenshaw et al.l 
I2OO0I : iDas et all |2005| ). the gas in the NLR is outfiowing 
along a hollow cone, with the SW side tilted toward us. Our 
line of sight is almost along, but outside, the near edge of the 
SW cone. As the cone seems to have a geometrical cutoff at 
~100pc from the nucleus, we are looking at the inner wall 
of the cone to the SW, while to the NE we are looking at the 
outer wall of the cone, which is at least partly behind the 
plane of the galaxy. We see higher excitation when we look 
at the inner wall of the cone and somewhat lower excitation 
when we look at the outer wall of the cone. As we do not see 
any reddening difference between the two sides of the cone, 
the only explanation seems to be that the ionizing radiation 
is attenuated more to the NE than to the SW. 

Recently, iKraemer et al.l (I2OO8I ) have mapped the ion- 
ization in the NLR of NGC4151 using the line ratio map 
[O lll]]/[0 11] obtained from the ratio of narrow-band HST 
images. They compare the observed emission-line ratios 
along the NLR with photoionization models and point out 
that, while the highest line ratios are observed along the bi- 
cone, there is also emission from regions outside the cone, 
with smaller line ratios, indicating lower excitation. They at- 
tribute this lower excitation to a weaker ionizing flux reach- 
ing these regions due to the presence of a "low-ionization 
absorber" which filter s the radiation. This a bsorber seems to 
be ionized gas, leading lKraemer et al.l (|2008l ) to propose that 
the low-ionization absorbers are dense knots of gas swept up 
by an accretion disk wind. 

Our excitation maps do not seem to show a decrease 
of the excitation perpen dicular to the bicone as observed 
bv iKraemer et al] ([200^). Instead, we find lower excitation 
to the NE than to the SW, as discussed ab ove. Maybe the 
differe nce between our results and those of iKraemer et ahl 
(|2008l ) can be attributed to the different optical depths 
probed by optical and near-IR observations. The difference 
in excitation between the NE and SW sides of the bicone 
could then be due to the presence of more absorbers be- 
tween the nuclear ionizing source and the emitting gas we 
see to the NE than to the SW. 



4.6 Physical conditions of the [Fell] emitting 
region 

There are many [Fe 11] emission lines in our near-IR spectra 
of the NLR of NGC4151, specially to the SW, where the 
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Figure 11. Model [Fell] emission-line ratios versus electron den- 
sity. Each sequence of models correspond to a different temper- 
ature - from bottom to top - 1, 3, 5, 10 and 20 X lO^K. The 
shaded regions correspond to the measured line ratio and their 
uncertainties within a region of aperture 0''3 X O'.'S centered in 
the [Fell] emission peak at 0'.'9SW from the nucleus. 



[Fe ll] emission is strongest. The near-infrared [Fe ll] emission 
lines are due to forbidden transitions between low energy 
levels, and their intensities are thus principally dependent on 
the electron density. We were able to measure four [Fe ll] line 
ratios in the region of strongest [Fe ll] emission (at « 0'.'9 SW 
of the nucleus) , from the fluxes listed in Table [T] which can 
be used as density indicators. We have performed a 16-level 
atom calculation for [Fe ll] and have obtained line intensities 
in order to compare with the observed ones in NGC 4151 and 
derive the electron density. Fig. l 111 shows the model results 
for four emission-line ratios as a function of the electron 
density for different values of the electronic temperature, 
along with the measured ratios and uncertainties. The line 
ratios indicate a density in the region of the strongest [Fe ll] 
emission of ~4000cm~"^. 

Most of the [Fe ll] emission-line ratios that are density 
indicators could be measured only at the location of max- 
imum [Fen] emission. The exception is [Fell]l. 533/1. 644, 
which could be measured at several other locations allow- 
ing the construction of an electron density map, shown in 
Fig |12l The uncertainty in the line flux of the 1.533/im line 
is quite large at the nucleus (see Table [!}, but decreases out- 
wards. It can be concluded from Fig. [12] that the electronic 
density in the [Fe ll] region decreases from the highest values 
at the nucleus lO^cm"^), to ~ 10^ cm"'^ within the in- 
ner f»0'.'5 from the nucleus. Outside this radius the density 
reaches the value of «4000 cm~^ obtained above. 

Deriving the electron temperature in the [Fe ll] emitting 
region is more difficult, as most observable emission lines 
arise from a limited range of upper level en ergies. Following 
iMouri. Kawara fc Taniguchil (|2000l ') and [McGregor et alj 
we looked for emission lines corresponding to 
transitions from the a*P term. There is a transition 
between the terms ''P4/3 and a.'^D-Y/2 corresponding to an 
emission line at 1.7489/im in the H-band. Although we 
have apparently detected this line, its uncertainty is too 
large to put a useful contraint on t he electron temperature. 
Another possibility suggested by iThompsoni (,1995. ) is to 




Figure 12. Distribution of electron density Ue obtained from the 
[Fell] 1.533/1.644 emission-line ratio, for an electron temperature 
in the range 5000K< Te <20000K. 



use the [Fe 11] 0.8619/im line strength, corresponding to the 
transition a'*P5/2^a'*F9/2 as an indicator of the electron 
temperature. Unfortunately, our observed spectral range 
does not cover this emission line and we had to look 
for previous observations a vailable in the literature in 
which this line was measured. lOsterbrock. Shaw fc Veilleuxl 
ljl990l ) have identified the [Fe 11] 0.8619/im emission 
line in their spectrum of NGC 4151, reporting a flux 
ratio [Fe 11] 0.8619/[S III] 0.9533=0.042. We have used 
this ratio and our measurement of the [S ill] 0.9533 
line flux at 0'.'9 SW from the nucleus in order to 
estimate the flux of the [Fell] 0.8619 emission line 
at this location, under the assumption that the 
[Fell] 0.861 9/[S III] 0.9533 ratio is the s a me a s that ob- 
tained by lOsterbrock. Shaw fc VeiUeuxl (|l990h . in spite 
of their much larger aperture. Using the fluxes from 
Tabled we obtain [Fe 11] 0.8619=2.09 x 10" ergscm"^s"\ 
and thus a ratio [Fe 11] 1.2570/0.8619=3.6. If we con- 
sider a reddening value E(B-V)~0.5, the corrected 
log([Fe 11] 1. 2570/0. 8619)=0.46±0.11 and the resulting 
temperature will be 15,000 ± 5000 K (considering the 
uncertainties in the line ratio combined with those in the 
reddening correction). 

The origin of [Fe 11] emission in galaxi es has been the 
subject of many studies, since the finding bv lForbes fc WardI 
d 19931 ) and lBlietz et"!!! l| 19941 ) of a correlation between the 
[Fell] and radio emission. This suggests that shock excita- 
tion by radio jets is a likely mechanism for produc tion of 
the [Fe 11] emission ()Dopita fc Sutherlandiri995l . ll996l '). How- 
ever, subsequent works have argued that the dominant exci- 
tation mechanism is photoionization, with shock excitation 
a ccounting for only ~20% of the [Fell] excitation in AGN 
(ISimpson et al.lll996l ). 

An important tracer of the origin of the [Fe 11] emission 
is the line ratio [Fe ll]/Pa/3, whose value ranges from <0.6 for 
starbursts to > 2 for supernova remnants, for which shocks 
are the dominant mechanism ( Ro driguez- Ardila ct alTI 
l2004h . Active galaxies have values for this ratio between 
0.6 and 2, suggesting that for ratios close to 0.6 photoion- 
ization is the dominant mechanis m, while for ratios close 
to 2, shock excitation dominates ([Storchi-Bergmann et al.l 



1 1999! : [Rodriguez- Ardila et al.1l2004f ). [Feii] /Pag ratios :^2 
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log [1(1,357 /^m), '1(0,862 ^m)] 

Figure 13. Model [Fell] line ratio 1.644/im/1.533/im (which 
depends mostly on the gas density) versus the line ratio 
1.257/im/0.862/im (which depends mostly on the temperature) as 
a function of density and temperature. The shaded region corre- 
sponds to the measured line ratios and their uncertainties within 
a region of aperture (/.'S X (X'S centered in the [Fe ll] emission peak 
at O'.'g SW from the nucleus. 



have been found by iRiffel et all (j2006bl ) in a near-IR IFU 
study the NLR of the Seyfert 2 galaxy ESO 428-G14 in re- 
gions showing a close association with emission knots of a 
radio image, supporting an origin for the [Fe ll] emission in 
shocks produced by the radio jet. 

In NGC4151, [Fell]/Pa/3 <1 within the inner C/.'6, in- 
creeising to values of ~2 outwards (see FigO, suggesting 
that photoionization dominates in the inner region while 
shocks dominate in the outer region. Indeed, inspection of 
the [Fell]/Pa/3 ratio map in Fig.[7]shows a relation with the 
radio structure: to the SW, the increase in the [Fe ll]/Pa/3 ra- 
tio does coincide with a flaring in the radio contours, which 
may have been produced by a shock between the jet and 
ambient gas, while to the NE, a similar increase is observed 
at the location of a radio knot. 

Another powerful indicator of the origin of the [Fe ll] 
em ission is the [Fenl 1.2570/[P ii] line ratio, as pointed out 
bv lOliva et all (|20o3). This line ratio is seldom used be- 
cause the [P ll] emission line is usually too faint. The qual- 
ity of our data has allowed the measurement of extended 
emission in this line. Besides having similar wavelengths, the 
[Fell] 1.2570 and [P ll] lines have similar excitation tempera- 
tures, and their parent ions have similar i onization potentials 
and radiative recombination coefficients. lOliva et all (|200ll ) 
have shown that, for a solar abundance, the above ratio is 
> 20, as observed in supernova remnants. But, in many as- 
tronomical objects this ratio is much smaller, because Fe 
is locked into grains. For example, in Orion, the above ra- 
tio is «2. In order to destroy the dust grains and release 
the Fe, fast shocks are necessary. Larger ratios than 2 indi- 
cate that shocks have passed through the gas destroying the 
dust grains, releasing the Fe an d thus enhancing its observed 
abundance. lOliva et al.l (j200l| ) have reported the first obser- 
vation of this line ratio in an extragalactic source, and here 
we report the first 2D map of this ratio in an extragalactic 
source. In NGC4151, [Fell]/[Pll] varies from 2 to 6, with 
the highest values observed in patches which do seem to be 



spatially correlated with the radio structure, as shown in 
Fig. [7] the highest ratios are observed at the location of a 
radio knot to the NE and at ~ O'.'S SW from the nucleus, 
where the flaring of the radio contours is observed. These lo- 
cations also coincide with those where the ratio [Fell] /Pa/3 
shows the highest values, supporting an increased contribu- 
tion of shocks to the excitation of [Fe ll] at these locations. 
Nevertheless, the [Fell]/[Pll] line ratio is never as high as 
«20, which suggests that the excitation mechanism is not 
only shocks, but includes photoionization. 



4.7 Physical conditions of the H2 emitting region 

As discussed by iRiffel et al.l l|2006bl ). the H2 emission 
lines can be excite d by two processes: (1) fluo rescence by 
soft-UV photons (|Black fc van Dishoeckl [iQSj) - present 
both in star-forming regions and around active nu- 
clei; or (2) thermal processes, produ ced either by X- 
ray (iMalonev. HoUenbach fc TielensI [199^1 or shock heating 
(|Hollenbach fc McKeelll989l i 

An emission-line ratio commonly used to investigate 
the origin of the H2 excitation is H2 2.1218/Br 7. In star- 
burst galaxies, where the main heating agent is U V radi- 
ation, H2/Br7 < 0.6 l|Rodriguez- ArdUa et all |2004| ) . while 
for Seyferts this ratio is larger because of additional thermal 
excitation by shocks or X-rays from the AGN. As observed 
in Fig.[7l this line ratio is <1 around the nucleus and along 
the axis of the bicone, but increases to values > 2 and up 
to 4 along the minor axis of the galaxy (and approximately 
perpendicular to the bicone axis). By comparing Fig[7]with 
Fig[6] which shows that the H2 intensity distribution avoids 
the region along the bicone axes, we propose that the low 
H2/Br7 ratio in these regions is due to the destruction of 
the H2 molecule by the strong radiation and particle fluxes, 
while the high values at the perpendicular direction could be 
due to fluorescence or to excitation by nuclear X-rays escap- 
ing in that direction or by shocks. In the case of NGC 4151, 
shocks in the region where H2 emission is observed, perpen- 
dicular to the radio jet, can be produced in gas flowing into 
the nuclear region or associated with star formation in the 
disk. It could be even possible that star formation is occur- 
ring in in-flowing gas. 

The H2 A2.2477/A2.1218 line ratio can be used to dis- 
tinguish bet ween th e rmal ( ~0. 1-0.2) and fluorescent (~0.55) 
excitation (iMouri 1993: iReunanen. Kotilainen fc Prietol 



l2002l : [Rodriguez- Ardila et al. II2004I ). The first line could be 
measured in a few locations along the NLR. The resulting 
line ratio does not vary much along the NLR, and the aver- 
age value is 0.13±0.02, supporting a thermal excitation for 
the H2. 

We can investigate further the excitation mechanism 
of H2 by using all its emission-line fiuxes in the K band 
to calculate its excitation temperature T^xc and the ratio 
of the ortho to para emission lines. As the H2 emission- 
line ratios seem not to show spatial variation, in order 
to improve the signal-to-noise ratio, we have measured 
the H2 fluxes within a O'.'S diameter aperture centered 
on the SE peak of the H2 intensity distribution. FoUow- 
ing lWilman. Edge fc Johnstone! (|2005l ) , we have investigated 
the relation: 
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log ( ' j = constant — (3) 

where Fi is the flux of the ith H2 hne, \i is its wavelength, Ai 
is the spontaneous emission coefflcient, gi is the statistical 
weight of the upper level of the transition and Ti is the 
energy of the level expressed as a temperature. This relation 
is valid for thermal excitation, under the assumption of an 
ortho : para abundance ratio of 3:1. T^xc (the reciprocal of 
the slope) will be the kinetic temperature if the H2 is in 
thermal equilibrium. 

The above relation is plotted in Fig. [TJ] where we have 
included the fluxes of most H[2 emission lines, excluding only 
the weakest ones (1-0 Q(6), 1-0 Q(7) and 2-1 S(l)). The re- 
sult is a tight straight line with T^xc =2155 K. The fact that 
all emission lines are on the relation confirms that the H2 is 
in thermal equilibrium at Texc, and that the rotational and 
vibrational temperatures are the same, ruling out a signifi- 
cant contribution from fluorescence. The fact that the ortho 
lines (1-0 S(l), 2-1 S(l), 1-0 Q(l), etc) give the same result 
as the para lines (1-0 S(0), Q(2), Q(4), etc) confirms that 
the ortho to para ratio is ~ 3 as assumed. This is also what 
is expected for thermal equilibrium. 

We have thus concluded that the H2 emission-line ratios 
are consistent with thermal excitation, which excludes UV 
fluorescence from the nucleus as a mechanism to excite the 
H2 emission. But there are two other possibilities: heating 
by nuclear X-rays or by shocks in the gas flowing to the 
nucleus. 

In order to test if X-rays from the AGN can account 
for the observed H2 line fluxes, we hav e used the models of 
iMalonev. HoUenbach fc Tielenj (119961 ) to estimate the H2 
flux emitted by a gas cloud illuminated by a source of hard 
X-rays with lu minosity I/jy. The ca lcula tions are performe d 
as described in lZuther et~all (|2007| ) and lRiffel (|2008l ). 
Using a X-ray flux of 4.51 x 10~^^ ergcm"^ s~^ and an ab- 
sorptio n column density of 7.5 x 10^^ cm^ for the nuclear 
source (jCappi et al.ll20()6l ). we obtain a X-ray flux at a typi- 
cal distance from the nucleus of 33 pc (0'.'5 from the nucleus, 
where H2 emission is observed) of 7.3 erg cm~^ s~^ and an 
effec tive ionization parameter of ^eff = 0.0 15. Using Fig. 6a 
from lMalonev. HoUenbach fc Tielei3 (11996'), we can then es- 
timate the resulting H2 flux at the Earth for an aperture 
corresponding to our pixel of O'.'l x 0'.'04 - which gives a 
solid angle of 9.6 x 10"^* sr - as 7 x lO"^'^ erg cm"^ s"\ 
Inspection of Fig[6] shows that the highest values observed 
are « 10~^® ergcm"^ s~^ and we conclude that excitation 
by X-rays emitted by the AGN can account for most of the 
observed H2 flux if a similar X-ray flux to that escaping in 
our direction also reaches the region of H2 emission. 

iPrestwich et al. I analysed near-infrared spectra 

of NGC 4151 obtained with > 5.4" diameter apertures over a 
ten year period from 1979 to 1989 and found that the II2 line 
emission was stable over this interval while the (broad) H I 
Br7 emission varied signiflcantly. This could argue against 
nuclear X-ray heating. However, the light crossing time of 
the H2-emitting region with its ~ 1" radial extent (Fig. [6]) 
is ~ 200 yr, so nuclear variations on decade timescales will 
be spatially diluted in large aperture measurements. Obser- 
vations of individual H2 clumps at higher angular resolution 
than possible with NIFS would be required to provide tighter 
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Figure 14. Relation between Nupp = and Eupp = Ti 

for the H2 emission lines for thermal excitation at temperature 
Texc =2155 K. Transitions plotted are from (left to right) 1-0 
Q(l), 1-0 S(0), 1-0 Q(2), 1-0 S(l), 1-0 Q(3), 1-0 S(2), 1-0 Q(4), 
1-0 Q(5), 2-1 S(l) and 2-1 S(3). 

constraints since the light crossing time of a 0.1"diameter 
clumps is still ~ 20 yr. 

What then is powering the H2 emission? Our observa- 
tions are consistent with two mechanisms: X-rays from the 
active nucleus and/or shocks in the inner galactic disk, pos- 
sibly prodAiced__bytlwaccretion flow along the large scale 
bar (iMundefl et al.lll995l ). which may be the origin of the 
H2 gas accumulated in the inner region. 

4.8 Mass of ionized and molecular gas 

The mass of ionized hydrogen can be estimated as Mhii = 
mp Ne Vhii, where A'^e is the electron density and Vhii is the 
volume of the emitting region. 

The product Vkii can be obtained from the expres- 
sion for the Br 7 flux obtained us ing the H I emission coeffi- 
cients listed in lOsterbrockl (|l989D : 

c o -7 w in-28 -^e ^HII -2 -1 , 

FBr7 = 2.7 X 10 — ^5 — erg cm s (4) 

where D is the distance to the galaxy in cm, units of density 
are cm~^ and we have assumed an electron temperature of 
10''K and density in the range 10^ < Ae < lO^cm"^. The 
resulting gas mass in solar masses is: 




The integrated Br7 fiux over the NLR is fa 4.2 x 
10~^* erg s~^ cm~^. As discussed above, there is an aver- 
age reddening of E(B-V)=0.5 al ong the NLR. Co rrecting for 
t his re ddening, using the law of lCardelli. Clayton & Mathi^ 
(| 19891 ') the Br7 flux is Fsv-y ~ 4.4 x 10"" erg s"^ cm~^. 
Adopting an electronic density of 100 cm^ we obtain Mhii ~ 
2.4 X 1O''M0. The adopted density is justified by the fact 
that the resulting H 11 mass, divided by the volume of two 
cones with height of 100 pc and opening angle of 75° - the 
approximate geometry of the observed intensity distribu- 
tions - does result in a density value of «100cm~^. It may 
be that the density is higher, and the filling factor is smaller 
than one, but the resulting mass will be the same. 

We can also calculate the mass of hot H2 (which emit 
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the K-band emissio n lines) as in iRiffel et al.l l|2008l ) and 
IScoville et al] l|l982h : 



H2 



2m,p FH2X2.121S 4-kD 



= 5.0776 X 10' 



V Mpc 



(6) 



where rup is the proton mass, Fff2A2.i2i8 is the line flux, 
D is the distance to the galaxy and is given in 

solar masses. For a typical vibrational temperature of 
Tuii, = 2000 K (similar to the value we have obtained) 
the population fraction is /i/=i,j=3 = f.22 x 10~^ and 
the transition probability is Aai-,-. = 3.47 x lO"'^ s~^ 



jRiffel et al. 20081: 


Turner, Kirbv-Docken & DalEarnoll 19771: 


IScoviUe et al. 1982 


). 



The total H2 2.1218 flux, integrated over the whole 
emitting region is 2.5 x 10~ ergs" cm~ . Correcting 
for the average reddening of E(B-V)=0.5, we obtain 
^"^2^2.1218 = 2.64 X 10~^* ergs~^ cm~^. Then, using the ex- 
pressions above we obtain Mjj2 =240 M0 . 

This value is 10* times smaller than that of H 11, but it 
should be noted that this mass is only of the hot H2, which 
emits because either shocks (probably from the accretion 
flow along the bar) or X-rays from the AGN excite the H2 
molecule. Most of the molecular gas in the nuclear region of 
galaxies is cold, wit h the hot-to-cold mass ratio ranging be- 
tween 10"^ to 10"^ (|Dale et al.ll2005h . Thus, the total mass 
(hot plus cold) of molecular gas is probably even larger than 
that of H 11. The presence of such a molecular gas reservoir 
around AGN is supported by radio observations of the cen- 
tral region of active galaxies, by, for example, the "NUGA" 
group, which report molecular gas masses in the range 10^- 
10^ Mq in th e inner few hundred parsecs of active galaxies 
jBoqne et al.ll2007l : iGarcia-Burillo et al. II2OO5I : iKrips et all 
20071. e.g.). In the case of NGC4151, as discussed above, 



Mundell fc Shon"3 (|l999l ) have measured inflows in radio ob- 



servations of H I along the large scale bar (at PA 130°), which 
may be the origin of a molecular gas reservoir accumulated 
around the nucleus, only a small part of which we observe, 
due to excitation of the H2 molecule. 



5 SUMMARY AND CONCLUSIONS 

We have mapped the emitting gas intensity distributions, 
reddening and excitation in the NLR of NGC4151 using 
Gemini NIFS observations of the inner 200x300 pc^ of 
the galaxy, covering the near-IR Z, J, H and K spectral 
bands, at a resolving power R^5000 and spatial resolution 
of ~8 pc. The main results of this paper are: 

• The intensity distributions in the recombination lines 
of H and He as weU as in the [S ill], [P 11], [ [Fe 11] and [S ill] 
emission lines, similarly to that in the optical [O ill] , are most 
extended along PA=60/240° - the axis of the previously 
known bicone. The emitting region is somewhat brighter and 
reaches a larger projected distance of «130pc from the nu- 
cleus to the SW (the near side of the bicone) than to the 
NE (the far side of the bicone) , where it reaches a distance 
of «100pc. The fluxes in the recombination lines decrease 
with distance r from the nucleus as oc , while those on 



the other emission lines above remain constant or even in- 
crease with distance from the nucleus along the NLR. 

• The H2 intensity distribution is completely different 
from that of the ionized gas, avoiding the region of the bi- 
cone, probably due to the destruction of the H2 molecule 
by the strong ionizing flux along the bicone. Most of the 
H2 emission seems to be coming not from the outflowing 
gas, but from gas which is in the plane of the galaxy. This 
is supported by its kinematics (Paper II). The concentra- 
tion of the H2 emission approximately perpendicular to the 
bicone indicates that there is attenuation of the ionizing 
flux at these locations, precluding the destruction of the H2 
molecule. This attenuation may be due to an obscuring torus 
and/or to the bottom pa r t of t he biconical outflow (as sug- 
gested bv lKraemer et al.l (|2008l )'). 

• The intensity distribution in the coronal line [Sivil] is 
resolved and that in the [Ca vill] is marginally resolved, con- 
sistent with an origin in the inner NLR. The intensity dis- 
tributions in these and the other coronal lines - [S vill] and 
[Six] - decrease steeply with distance as oc , similarly 
to those of the calibration stars. 

• The line ratios [Fe 11] Al. 257/1. 644 and Pa/3/Br7 were 
used to map the reddening along the NLR, which ranges 
from E(B-V)=0 to E(B-V)=0.5. Within O'.'S from the nu- 
cleus, the reddening increases up to E(B-V)S^1. 

• The [Slll]/Pa/3 line ratio has a value of ~11 to the 
SW (typical for the NLR of Seyfert galaxies), where we are 
looking at the inner wall of the near cone, and is « 30% lower 
to the NE, where we are looking at the outer wall of the far 
cone. This difference seems not to be due to reddening and 
indicates lower excitation to the NE. 

• The line-ratio map [Fe 11] 1.257/[P 11] 1.187 of the 
NGC4151 NLR is the first 2D such map of an extragalactic 
source, and, similarly to the [Fe 11] 1.257/Pa/3 ratio shows 
larger values at the locations where the [Fell] emission is 
enhanced at « 1" from the nucleus. The increase in these 
line ratios maps the shocks produced by the radio jet in 
the NLR, which release the Fell usually tied up in dust 
grains. This is confirmed by the correlation between the 
linhe-ratio maps and the outer parts of the radio map. From 
the many emission lines observed at the locations where the 
[Fell] emission is enhanced, we have obtained the gas den- 
sity, Ne «4000cm"^ and temperature T^ « 15000±5000K. 

• From the fluxes of 10 H2 emission lines, we have con- 
cluded that the H2 emitting gas is in thermal equilibrium 
at the excitation temperature Tf^xc = 2155 K, ruling out 
any significant contribution from fluorescence to the exci- 
tation of the H2 molecule. The thermal excitation may be 
due to X-rays from the AGN escaping perpendicular to the 
bicone axis or else to sh ocks produced by the accretion flow 
observed along the bar (|Mundell et al.|[l995l ) when the gas 
reaches the nuclear region. 

• We have calculated the mass of the ionized and molec- 
ular gas, obtaining for the former Mhii ~ 2.4 x 10^ Mq 
and for the latter only ~ 240 M© . This small mass is 
nevertheless only that of the "hot skin" of a probably much 
larger molecular gas mass. 

• The distinct intensity distribution and smaller temper- 
ature (as well as distinct kinematics; see Paper II) of the H2 
emission when compared with those of the ionized gas, sup- 
ports a distinct origin for the emitting gas. The H2 emission 
is probably a tracer of a large molecular gas reservoir, built 
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up from the accretion flow observed along the large scale 
bar, which may be the source of fuel to the AGN. The H2 
emission can thus be considered a tracer of the feeding of 
the AGN, while the ionized gas emission, which maps the 
outflowing gas, is a tracer of the feedback from the AGN. 
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